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ABSTRACT 
The  stresses  near  the  root  of  an  F8U-3  airplane  wing  loaded  in 
torsion  were  found  using  the  matrix-force  method  of  analysis.   The  wing 
was  represented  by  a  simplified,  idealized  model  consisting  of  axial- 
load  carrying  flange  bars  and  shear  carrying  cover  plates.   The  theoret- 
ical stresses  in  the  idealized  model  were  found  by  using  the  IBM  360/67 
computer  to  solve  for  the  redundant  loads.   These  stresses  show  favor- 
able comparison  with  the  experimental  stresses  determined  from  strain 
gages  mounted  on  the  cover  skin  of  the  wing.   This  study  shows  that  the 
matrix-force  method  is  applicable  to  thick-skinned  wings. 
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INTRODUCTION 

The  purpose  of  this  study  is  to  predict  the  stresses  in  an  F8U-3 
airplane  wing  loaded  in  torsion  using  the  matrix-force  method  of  analysis. 
The  stresses  determined  analytically  are  compared  with  the  stresses  ob- 
tained experimentally  in  order  to  determine  the  accuracy  and  consistency 
of  the  results. 

The  primary  structure  of  the  F8U-3  wing  center  section  consists  of 
a  simple  multicell  torque  box,  as  shown  in  Figure  1.   Spanwise  bending 
is  carried  by  the  relatively  thick  skins  which  are  fabricated  from 
7079-T6  aluminum  alloy  plate.   The  skins  are  machined  to  varying  thick- 
nesses both  spanwise  and  chordwise  in  order  to  obtain  an  optimum  distri- 
bution of  material.   The  skins  are  continuous  from  the  wing  fold  rib 
inboard  to  the  wing  centerline  and  from  the  front  beam  to  the  rear  beam. 
It  is  this  area,  between  wing  fold  ribs,  port  and  starboard,  and  Front 
and  Rear  Beams,  fore  and  aft,  which  defines  the  wing  center  section, 
the  area  of  interest  in  this  study.   The  sweep  of  the  wing  is  42°  at 
the  quarter  chord.   The  wing  thickness  varies  from  47«  at  the  tip  to 
5.017o  at  the  root.   It  is  this  combination  of  a  high  degree  of  sweep 
and  low  thickness  ratio  which  raises  problems  with  structural  analysis 
of  the  wing.   Due  to  bending-torsion  interaction,  or  sweep-coupling, 
structural  analysis  by  engineering  beam  theory  methods  runs  into  diffi- 
culty when  analyzing  stresses  near  the  wing  root. 

This  thesis  describes  the  matrix-force  method  of  analysis  for 
determining  the  stress  distribution  in  such  a  wing.   Among  the  advan- 
tages of  this  method  are: 

(1)   It  results  in  a  consistent  set  of  internal  stresses  and  external 
deflections . 
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Actual    wing 


WS  81.98 


Idealized  wing 


FIG.  I 
F8U-3     WING    CENTER     SECTION 


(2)  The  structure  can  be  altered  and  re-analyzed  with  only  minor 
changes  in  the  computer  data  cards. 

(3)  The  "Bookkeeping"  can  be  simplified  and  organized  in  a  logical 
manner . 

The  matrix-force  method  of  analysis  uses  forces  as  the  unknowns  and 
treats  the  structure  as  an  assembly  of  elastic  components.   The  method, 
as  described  in  Reference  1,  is  divided  into  three  parts: 

A.  Idealize  the  structure 

1.  Replace  the  actual  structure  with  a  mathematical  model. 

2.  Make  the  model  statically  determinate  by  cutting  redundant 
members . 

3.  Write  the  equilibrium  equations  that  determine  the  load  distribu- 
tion in  the  structure 

B.  Calculate  and  tabulate  the  flexibilities  of  the  individual  members 
comprising  the  structure. 

C.  Perform  matrix  operations. 

In  this  report,  the  wing  is  "cut"  outboard  of  the  intermediate 

rib  (Y  =  81.98),  and  it  is  assumed  that  the  stress  distribution  here 
w 

can  be  accurately  calculated  using  the  engineering  theory  of  torsion, 
T  —  ^2.Aq .  The  dihedral  of  the  wing  is  neglected,  and  the  wing  is 
assumed  to  be  supported  on  four  rigid  supports  located  at  A,  B,  C, 
and  D  in  Figure  2.   The  aft  supports  are  pins  which  can  support  fore 
and  aft  and  vertical  loads.   The  forward  supports  can  support  vertical 
loads  only.   The  only  loading  condition  considered  in  this  study  is 
symmetrical  wing  torsion.   From  symmetry,  the  reactions  at  all  four 
support  points  are  determined  by  equilibrium  as  functions  of  the  applied 
loads . 
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The  model  selected  to  represent  the  wing  structure  consisted  of 
flat  plates  which  carry  only  shear  and  cylindrical  bars  which  carry 
only  axial  load.   This  model  was  selected  for  compatibility  with  the 
previous  work  done  in  Reference  2.   Analysis  of  the  model  can  be  as 
accurate  as  the  matrix-force  method  and  arithmetic  permit;  however, 
the  correlation  between  theoretical  and  experimental  stresses  will 
depend  upon  how  closely  the  model  represents  the  actual  structure. 
This  study  is  not  an  attempt  to  define  the  best  model  for  the  analysis, 
but  to  calculate  stresses  in  what  was  tacitly  assumed  to  be  an  adequate 
model.   Most  of  the  examples  of  similar  studies  currently  available  in 

the  literature  deal  with  thin-skinned  shells.   There  is  little  evidence 
on  whether  or  not  such  models  are  suitable  for  the  thick-skinned  wing 
plates  used  in  this  study.   An  important  factor  in  constructing  the 
idealized  model  is  the  change  of  the  wing  geometry  to  model  geometry. 
The  model  consists  of  bars  and  flat  shear  carrying  plates.   It  is  ex- 
tremely important  to  insure  that  these  plates  are,  in  fact,  mathemat- 
ically flat.   For  example,  the  depths  of  the  ribs  at  WS  25.095  and  WS 
81.98  were  such  that  the  four  points  AFGH  in  Figure  1  did  not  lie  in  a 
plane.   Therefore, panel  AFGH  is  warped  and  has  panel  edge  shears  which 
have  components  perpendicular  to  plane  AFGH.   These  shears  are  not 
accounted  for  in  the  equilibrium  equations  written  for  a  flat  panel 
and  will  produce  errors  in  the  solution.   Reference  3  shows  that  for 
a  similar  wing  such  discrepancies  yielded  determinate  member  loads 
which  were  out  of  equilibrium  as  much  as  5  to  10%.   After  altering 
the  beam  depths  to  make  the  panels  flat,  this  trouble  disappeared.   The 
alteration  required  a  change  of  less  than  1%  of  each  dimension. 
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After  determining  model  geometry  from  the  adjusted  dimensions  of 
the  wing,  the  cross  sectional  area  of  the  bars  and  the  thickness  of  the 
shear  panels  were  calculated.   Each  bar  area  of  the  model  consists  of 
the  area  of  the  beam  flange  plus  one  half  the  skin  area  between  adjacent 
beams.   The  beam  web  areas  were  not  included,  using  the  methods  of 
Reference  4.   The  effect  of  combining  the  axial  load  carrying  skin  with 
the  beam  flange  to  form  the  model's  cylindrical  bar  is  to  remove  Poisson's 
effect  from  the  analysis.   This  means  that  any  chordwise  contraction  of 
the  cover  skins  is  neglected.   The  wing  skin  thickness  varies  in  the 
chordwise  direction  in  section  ABCD  and  in  both  spanwise  and  chordwise 
directions  in  section  AJBEF.   An  average  skin  thickness  between  beams 
was  determined  and  used  as  the  shear  plate  thickness.   Using  constant 
thickness  shear  plates  in  place  of  tapered  skin  simplifies  the  analysis 
but  introduces  possible  sources  of  error.   Model  geometry  and  areas  are 
listed  in  Tables  I  and  II.   Diagrams  of  cross  sectional  geometry  are  shown 
in  Figure  3,  4,  and  5.   An  exploded  view  of  the  outer  section  of  the 
idealized  model  from  Y  =  25.095  to  Y  =81.98  is  shown  in  Figure  6. 

WW 

An  exploded  view  of  the  intermediate  rib  at  Y  =81.98  is  shown  in 

w 

Figure  7. 

2.   Make  the  model  statically  determinate  by  cutting  redundant  members, 
The  most  complex  beam  that  is  truly  determinate,  so  that  the  flange 

and  shear  loads  can  be  determined  from  equilibrium  equations,  is  the 

3- flange,  3  web  beam  shown  in  Figure  8. 
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Fig  8. 

DETERMINATE  BEAM 
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In  this  beam,  loads  can  be  obtained  from  equilibrium  only,  without  con- 
sidering flange  areas,  material  stiffness,  or  deformation  of  the  structure. 
In  order  to  make  the  idealized  model  of  this  study  a  determinate  structure, 
it  was  necessary  to  cut  the  5  interior  webs  and  11  of  the  flanges,  re- 
sulting in  the  structure  shown  in  Figure  9. 


■Flange    14 


Fig.  9 
CROSS  SECTION  OF  DETERMINATE  SYSTEM 

The  number  of  redundancies  arising  from  the  geometry  of  a  box  structure 
has  been  given  by  Reference  4  and  others  as  B  +  N  -  4  where  B  is  the 
number  of  longitudinal  effective  flanges  continuous  across  each  junction 
and  N  is  the  number  of  closed  cells  stiffened  by  ribs  at  tip  end  of  bay. 
For  the  outer  section  of  this  model,  B  =  14   and  N  =  6, which  results  in 
16  degrees  of  redundancy.   This  agrees  with  the  16  members  which  were 
cut  to  provide  the  determinate  structure  of  Figure  7. 

3.   Write  the  equilibrium  equations  that  determine  the  load  distri- 
bution in  the  structure. 

Considering  the  outer  wing  section  from  station  25.095  to  81.98, 

there  are  the  following  49  unknown  member  loads. 

12  skin  panels 

14  beam  flanges 

7  beam  webs 

10  rib  flange  loads  at  station  81.98 

6  rib  webs  at  station  81.98 

49    total  unknowns 
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The  redundant  forces  come  from  the  following  members 


1 

F.B.    , 

,   upper  flange 

2 

FR.I.B. 

,   both  flanges 

2 

F.I.B. 

both  flanges 

2 

C.I.B. 

both  flanges 

2 

A.I.B. 

both  flanges 

2 

R.I.B. 

both  flanges 

5 

webs  8  • 

•  12 

16    redundant  forces 

Thus,  there  are  49  -  16  or  33  equations  of  equilibrium  required  to 
characterize  the  model.   These  equations  are: 

a)  Six  equations  of  equilibrium  of  free  body  of  bay  between  stations 

Y  =  25.095  and  81.98 
w 

(1)  P  =  0 

x 

(2)  P  =  0 

y 

(3)  p  =  o 

z 

(4)  M   =0 

xx 

(5)  M   =0 

yy 

(6)  M   =0 

zz 

b)  Eleven  equations  of  equilibrium  of  the  redundant  beam  flanges,  elements 

16,  18,  20,  22,  24,  29,  31,  33,  35,  37,  and  39. 

c)  Ten  equations  of  forces  and  moments,  F  =  0  ,     M   =0   of  the 

x        yy 

intermediate  beam  segments  of  the  rib  at  Y  =81.98 

w 

d)  Six  equations  of  equilibrium  of  the  rib  posts  on  the  beam  ends  at 

Y  =81.98,  elements  53  through  58. 

w 

The  addition  of  the  center  section  Y  =  0  to  25.095  can  be  thought 

w 

of  as  an  extention  of  the  horizontal  flanges.   The  forces  in  the  X  and  Z 

directions  are  reacted  by  the  horizontal  and  vertical  supports  at  Y  = 

25.095,  but  the  Y  forces  carry  through  to  the  centerline,  Y  =  0. 

w 
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Although  the  center  section  is  16°  redundant  (the  same  as  the  out- 
board section) ,  these  are  not  new  redundancies  but  meraly  components 
of  the  forces  in  the  outer  section.   Further  investigation  shows  that  for 
symmetrical  wing  loading,  the  bending  meoment,  M  ,  is  constant  from 

Y  =  25.095  to  Y  =  25.095  and  that  torsional  and  vertical  shear  are 
w  w 

zero.   The  axial  load  on  each  flange  at  Y  =  25.095  is  the  same  as  the 

°      w 

load  on  each  flange  at  Y  =  25.095.   Hence  there  is  no  shear  on  any 

w 

panel  or  web  in  the  center  section.   (It  was  later  shown  that  due  to 
a  slightly  unsymmetrical  load,  there  are,  in  fact,  shearing  stresses 
in  the  center  section  panels.)   Since  the  addition  of  the  center  section 
results  in  no  new  redundancies  and  no  new  external  loads,  there  can  be 
no  more  than  the  33  independent  equilibrium  equations  already  proposed 
for  the  outer  section.   Since  there  are  more  than  33  structural  elements 
involved,  there  are  several  possible  combinations  of  33  independent 
equations.   Those  chosen  were  considered  to  be  those  most  likely  to 
have  the  larger  stress  concentrations. 

B.   Calculate  and  tabulate  the  flexibilities  of  the  individual  members 
comprising  the  structure. 

As  described  in  Section  C,  the  internal  strain  energy  of  the  wing 
may  be  written  as: 

Z£U  -  L«.  a.  •■■  o%  J  [%,4  J  «1 1  (1) 

The  elements  in  the  F  matrix  define  the  elasticity  of  the  structure, 
If  skin  thicknesses  or  flange  areas  are  changed,  only  the  elements  of 
the  F  matrix  will  be  affected.   If  the  geometry  of  the  model  is  changed, 
the  equations  of  equilibrium  will  also  be  affected.   Multiplying  out 
equation  1, 
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2EU  =    *.f„  +  <*z  /L  *   ••  a,  Fhn  +2q,qx  Fiz   1-  2a,a0  Fl3  + 


(2) 


Each  term  of  Equation  2  can  be  regarded  as  a  small  increment  of  strain 

2 
energy  ^  U     associated  with  the  load  a.   or  a.  .  a..   Since  the  load  a. 

l      i    j  l 

may  be  applied  to  more  than  one  structural  element,  the  corresponding 

F..  term  may  be  made  up  of  contributions  from  several  structural  elements 
n 

The  F  matrix  must  be  a  square  symmetric  matrix  where  F.  .  =  F... 

1.   Flange  Bars. 

The  calculation  of  the  F  matrix  elements  for  the  flanges  was  done 
following  the  example  of  Reference  4.   The  energy  in  a  uniform  bar  under 
a  linearly  varying  axial  force  is 


u  = 


/ 


2AE 


dx 


(3) 


where    o     and   q -  are  end  forces   as   shown  in  Figure   10. 
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?- 


D -R 


£ 


Fig.  10 

FIANGE  BAR 
Rewriting  Equation  2  for  this  single  element 

2  *  2 


2EU  =%,£,+ ?%,  %*  fm    +  %z   £ 


(4) 


therefore 


f"        L   ho3; 


2L 
L 


dx 


(5) 


T"       3A  ** 


(6) 
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"M. 


-ML 

and 


(7) 
f  =   f    = 


''*  "  'zi        6A 

In  the  case  of  linearly  tapered  members,  the  flexibilities  may  be 
obtained  as  functions  of  the  end  area  ratios.   Reference  1  shows  the 
flexibility  influence  coefficients  modified  by  a  function,  Q,  of  the 
end  areas  which  was  obtained  from  a  graph  of  Q  versus  A.  /A.,  the  end 
area  ratio. 

£  .  -i-  d>  f  =—  <&  f  =  -^  6  (8) 

For  most  of  the  flange  elements,  the  generalized  end  loads,  q. 
and  q.  were  linear  functions  of  the  unkown  member  loads,  or: 

?.  =  <«,  *<*,  ^  "  •<"4i,  (9) 


?j  "  0  *'  * .ft  «e  '  '  '  '  ft  ^  (10) 

where  C&*    =  the  coefficient  of  the  i'th  member  load  applied  to  one  end 

of  the  bar, 

and  Q.    =  the  coefficient  of  the  i'th  member  load  applied  to  the  other 

end  of  the  bar. 
Now  taking  each  force  acting  on  the  ends  of  the  flange  individually, 
that  is: 

f ,  =  °t>  a<  (ii) 

?*   *    P>  *>  (12) 

and  substituting  into  Equation  4,  we  find  that 

beu  =  a#Y<£  +eu.p,fa  y  (3/4  )  d3) 
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and  from  Equation  1,  for  a  single  force,  a 

2EU  =    a?  F„ 

Comparing  Equations  13  and  14, 

Similarly,  for  the  other  forces  acting  on  the  same  flange, 


22. 


and  for  the  cross  product  term, 


(14) 


(15) 


(16) 


(17) 


The  forces  on  element  14  are  shown  in  Figure  11. 

For  example,  the  flexibility  influence  coefficients  for  element  14  were 

determined  as  follows: 


%  * 


-! 

►- 

A,  =  /.  3 

-t                     -.._/—    TO    "?                                   «_ 

Az--/.3 

—  i.      /y.  j                  — ■» 

Fig.    11 

FORCE   DIAGRAM  OF  FLANGE  BAR 

*»  -  h,  A  = *«  ^  -  ^ 

y 

«=Z7  &.  =  /0- /67 

?/  =  -3.23/ a^  +ax/-/o.3S3a. 


?*-     -^e  **>**„ 


2  2, 

/\      =*0.  233 (-/0. 9*3)  *  *  *  /o./il^/o.  5si)(o)  +ZO,  333 (o)    =  2439.  433 
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For  element  14,  there  are  4  more  coefficients  to  be  found;  F 


F21,  21'  F2,  21'  and  F20,  21'   The  F2,  2 


20,  20' 
coefficient  for  element  14  is 


added  to  other  F   ~  coefficients  for  elements  7,  39,  101,  and  126  to 
make  up  the  total  F9   _  flexibility  influence  coefficent  found  in  the  F 
matrix. 

2.   Symmetrical  Shear  Panels 

Reference  3  gives  the  energy  of  the  symmetrical  trapezoidal 
panel  in  terms  of  the  shear  force  on  the  shorter  parallel  side,  as: 


2E 
so  that 


vH^^lH^'W 


a 


may  be  written  for  F   -  of  element  1 


'-'    (w. 


33S(/2,46o) 


S/£  (.2oo) 


■g  (.OOG78)    +  /32 


The  forces  on  element  1  are  shown  in  Figure  12 


W, 


a  =.8Z7a, 
c  1 


Fig.  12 
FORCE  DIAGRAM  OF  BEAM  WEB 


(18) 


(19) 


Other  F1   1  coefficients  appear  from  elements  40  and  51  and  add  to  pro- 
*■  >  *■ 

duce  the  final  F     found  in  the  F  Matrix. 
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The  shear  panel  edge  loads  were  determined  by  assigning  a  unit 
shear  force,  a.,  to  the  outboard  edge  of  all  panels  and  the  forward 
edge  of  the  rib  webs.   From  the  equilibrium  of  the  free  body,  forces 
a,  ,  a  ,  and  a  may  be  found  as  functions  of  the  unit  force  a.. 

%  -   ai  h   '   Wc 

a  =  a  W.  /  W 
c    11    c 

a,  =  a.  1   /  W 
d     id    c 


(20) 


The  geometry  of  the  panels  and  the  load  ratios  are  shown  in  Table  III. 

3.   Non-symmetrical  shear  panels. 
a)   The  cover  skin  panels  have  roughly  the  shape  of  a  parallelogram, 
with  the  inclined  sides  skewed  slightly.   The  maximum  difference  in 
inclination,  0-8=   3°  34',  occurs  at  the  forward  panels,  elements  38 
and  15.   These  panels  were  treated  in  accordance  with  Reference  6 
which  shows  that  the  energy  of  the  swept  panel  is  given  as 


££CT  = 


r/^-^K^37fe))] 


where   T  =  Cot  ~@  +  cot  8  cot  @  +  cot  * 6 

The  forces  on  element  15  are  shown  in  Figure  13. 


(21) 
(22) 
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so  that 

'  nsr\ 
may  be  written  for  F 


Fig.    13 
FORCE  TJIAGRAM  OF   COVER   PANEL 


=    E 


2 

Za4z 


20,    20 


of   element    15. 
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b)   The  cover  skin  panels  introduce  additional  forces  due  to  interaction 
between  torsion  and  bending.   Twist  of  a  swept  wing  produces  bending  and 
vice-versa.   This  interaction  was  introduced  by  considering  the  stress 
distribution  in  the  rhomboidal  cover  skin  panels,  as  shown  in  Figure  12. 
These  panels  were  assumed  to  be  parallelograms  in  order  to  simplify  the 
analysis.   Consider  the  plate  shown,  loaded  with  a  uniform  shear  force, 


a  ,  and  axial  load,  f  . 


a 


-  a 


f  =  2a   tan  A 
Fig.  14 
SWEEP  COUPLING  OF  COVER  PANEL 
The  axial  force,  f  ,  was  assumed  concentrated  at  the  flanges  which  were 
made  up  mostly  of  skin  material,  so  that  the  cover  panels  carry  only 
shear.   However,  it  is  seen  in  Figure  12  that  the  shear  forces  along 
the  sides  of  the  panel  produce  additional  tensile  stresses,  f  ',  which 
do  work  on  the  displacements  produced  by  the  primary  axial  stress,  f  . 
For  this  panel,  the  strain  energy  can  be  written. 


2  2 


(23) 


The  third  term  in  Equation  22  is  due  to  the  interaction  between  the 
axial  and  shearing  forces  and  does  not  appear  in  an  equation,  written 
for  a  rectangular  panel.   This  term  causes  the  stress  distribution  in 
this  model  to  show  a  root  stress  concentration.   Reference  3  states  that 
"Studies  made  on  simple  one  cell  swept  box  models  with  streamwise  ribs, 
utilizing  parallelogram  shaped  cover  sheets,  show  that  unless  this  cross- 
product  term  is  introduced,  there  will  be  nothing  in  the  analysis  to 
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differentiate  between  a  straight  and  swept  wing,  and  there  will  be  no 

3 
root  stress  concentration  factor"  .   For  the  purpose  of  calculating  the 

energies  in  the  flanges  adjacent  to  each  panel,  this  axial  load,  f ' 

t  > 

was  aplit  equally  and  added  to  the  end  of  each  flange.   These  coupling 

loads  are  shown  as  dashed  vectors  in  Figure  6. 

C.   Perform  matrix  operations. 

The  matrix-force  method,  as  described  in  References  1  and  3, 

expresses  the  energy  of  the  entire  structure  as  the  successive  product 

of  three  matrices. 

3£U   -  \a\  [f\  {a}  <2« 

where  a,  a  are  the  internal  loads  acting  on  the  members.   The 

In 

F  matrix  is  symmetrical  and  represents  the  elasticity  of  the  structure. 
The  internal  loads  were  found  by  writing  the  equilibrium  equations  for 
the  determinate  structure  and  for  each  redundant  member.   The  number 
of  equations  must  equal  the  total  number  of  unknown  internal  loads  less 
the  number  of  redundants.   The  equations  are  written  so  that  all  the 
independent  external  loads  and  redundants  are  collected  on  one  side  and 
non-redundant  loads  on  the  other  side  of  the  equality; 
such  as 

M  («*]  =  [£]     {P\ar}  HP 

where  A  and  B   are  matrices  of  the  coefficients  of  the  loads. 

Solving  for  the  non-redundant  loads,  a    ,  as  functions  of  the  applied 

nr 

loads,  P  ,  and  the  redundants,  a   ,  we  get: 

{anr}  =  ^[B]{P\ar}  =  [j}{p\ar}  9» 


33 


Then,  by  adding  the  trivial  relation  between  the  redundants 


and  expanding  J  to  include  Equation  4  results  in 
which  may  be  transposed,  yielding  Equation  6 

{«}  -  l«j  -  k  i  aA  Wr 

Now,  substituting  into  Equation  1, 

2EU=   [p\ar\6TFG    {P\CLr} 


(27) 


(28) 


(29) 


(30) 


H 


is  a  symmetric  matrix.   Its  members  are  the  coefficients  of  the 


strain  energy  terms  when  the  energy  is  written  as  a  function  of  applied 


and  redundant  loads.   Partitioning 


H 


into   four   sub-matrices, 


[«] 


"pup  '£    pxr 


H. 


2/ 


r    | 


H 


22      rxr 


yP  -  appl /ed   loads  -  33 
r  -    redundant  /oads  -  /6 


Rewriting  Equation  29. 


2£V*    P  H„   +2Par  Hl£+ar  H& 


(31) 


Minimizing  the  strain  energy  with  respect  to  the  redundant,  a  ,  yields 


12  r     2Z 


Rearranging 


M-  - 


en]  {p} 


(32) 


(33) 
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The  matrix  of  all  the  unknown  internal  loads  is  given  by 


total 


(a) 

m     - 

TOTAL 


P*f> 


£A/. 


W 


r  J 


n 


(34) 


S   is  the  unit  load  distribution  matrix  for  the  indeterminate  structure. 

Each  S.   term  is  the  value  of  the  i'th  member  load  when  the  structure 

1m 

is  subjected  to  a  unit  value  of  the  m'th  applied  load. 

Using  the  notation  of  Reference  5,  the  strain  energy  can  be 
written: 


U  =    ±[f\[c]{p] 


(35) 


where  C   is  the  matrix  of  flexibility  influence  coefficients, 
stituting  from  Equation  11  into  Equation  23, 


Sub. 


2E 


P 


J[#]MH 


therefore 


C 


■it 


rr 


F 


-     / 


I  \EN 


GTG 


=   _L     EN   H    EN 


and  substituting  EN  =  — 


H 


1-1 


22 


H 


21 


from  Equation  10, 


c  =  ^lk 


H 


IZ 


H. 


12. 


-I 


H 


Note  that  Equation  34  may  be  expanded  to  read 


a.  =  s.,  P.  +s..  P,  ...  s.   P 
1     111    122      1mm 


a2  =  S21  Pl  +  s22  P2 


s9   P 
2m  m 


(36) 


(37) 


(38) 


(39) 
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Therefore,  each  internal  load,  a.  ,  may  be  found  by  substituting  a 
particular  set  of  applied  loads,  ^Pj  .   For  any  other  loading  condition, 
only  the  matrix  of  applied  loads  needs  to  be  changed  to  obtain  the  new 
internal  loads  for  that  loading  (flight)  condition.   This  enables  a 
rapid  evaluation  of  the  maximum  load  on  any  element  over  a  range  of 
applied  loads. 

A  computer  program,  originally  set  up  by  Reference  2  for  the  CDC  1604 
computer  was  rewritten  in  Fortran  IV  to  fit  the  IBM  360/67  computer  and 
appears  in  Appendix  A.   This  program  performs  the  complete  matrix  oper- 
ation outlined  in  Section  C  with  input  A  and  B  from  the  flexibility  matrix 
and  the  applied  loads,  P.   Matrices  A,  B,  F,  and  P  appear  in  Tables  IV, 
V,  VI,  and  VII. 
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EXPERIMENTAL  RESULTS 

In  order  to  ascertain  whether  or  not  the  results  obtained  from  the 
strain  energy  method  were  accurate  for  this  type  of  thick  skinned  wing 
construction,  it  was  necessary  to  have  experimental  data  with  which  to 
compare  the  results  of  the  idealized  model.   The  wing  is  mounted  in  a 
jig  which  closely  approximates  the  actual  method  of  joining  the  wing  to 
the  aircraft  fuselage,  as  shown  in  Figure  15.   Since  the  wing  is  a 
continuous  structure,  the  root  fixity  coefficient  is  not  changed  as  when 
a  single  span  wing  is  mounted  on  a  completely  rigid  bulkhead.   The 
loading  condition  investigated  was  a  torque  of  336,000  in.  lbs.   This 
load  was  considered  to  be  well  within  the  elastic  range,  but  still  of 
sufficient  magnitude  to  give  adequate  strain  levels  throughout  the 
structure  for  repeatable  strain  measurements.   Some  strain  gages  were 
installed  on  the  wing  when  received  from  Ling-Temco-Vought  Aerospace 
Corp.,  and  others  were  installed  at  the  Naval  Postgraduate  School. 
Figure  16  shows  the  wing  with  one  cover  skin  removed  for  installation 
of  strain  gages.   Wing  section  Y  =  34.5  was  particularly  well  instru- 
mented, with  strain  gage  rosettes  placed  on  both  sides  of  top  and  bottom 
skins  to  check  for  differential  bending  effects.   Strain  gage  data  were 
reduced  to  principal  stresses  and  directions  through  computer  program 
"Recrose",  written  in  Fortran  IV  for  the  IBM  360/67  computer,  as  shown 
in  Appendix  B.   The  torque  loads  were  applied  with  four  hydraulic  cylin- 
ders through  a  heavy  fitting  attached  to  the  wing  fold  lugs.   The  applied 
torque  was  designed  to  be  perpendicular  to  the  elastic  axis.   The  load 
was  applied  at  both  wing  fold  axes;  and  to  ensure  a  symmetrical  load, 
the  four  hydraulic  cylinders  were  supplied  through  a  common  manifold. 
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Load  magnitude  was  determined  by  a  manifold  pressure  gage  in  the 
hydraulic  pressure  system  and  checked  by  four  Dillon  dynamometers 
attached  to  each  cylinder.   The  strain  readings  were  taken  with  a 
Budd/Datran  digital  strain  indicator  (Model  TC22) .   All  gages  were 
wired  to  a  common  junction  panel  as  shown  in  Figure  17.   All  gages 
were  balanced  at  a  zero  reading  before  applying  the  load.   Eighteen 
gages  were  read  at  each  loading  in  about  one  minute.   This  minimized 
the  chance  for  drift  in  gage  reading.   All  gages  were  rechecked  for 
their  zero  reading  after  the  load  was  released  and  sufficient  time 
had  elapsed  for  hysteresis  effects  to  be  minimized.   All  testing  was 
done  at  night  to  eliminate  the  thermal  effects  of  sunlight  shining  on 
the  wing.   Reference  2  contains  all  strain  gage  coordinates,  gage  fac- 
tors and  resistances.   It  also  shows  that  the  strains  are  in  the  linear 
region  by  comparing  the  gage  readings  at  a  loading  of  8000  psi  (336,000 
in.  lbs.  of  torque)  with  those  taken  at  7000  psi  (294,000  in.  lbs.  of 
torque)  and  extrapolated  to  8000  psi. 

Strain  readings  for  the  rosettes  on  the  cover  panels  appear  in 
Table  VIII.   The  readings  were  averaged  for  the  back- to-back  gages 
and  entered  in  program  "Recrose"  to  compute  the  stresses.   The  normal 
and  shearing  stresses  and  the  principal  axes  for  the  skin  panels  are 
also  shown  in  Table  VIII. 
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DISCUSSION  OF  RESULTS 

There  are  two  areas  of  comparison  in  this  study.   The  first  is  the 
comparison  of  the  theoretical  results  with  the  experimental  results,  and 
the  second  is  the  comparison  of  the  theoretical  analysis  of  the  wing 
with  and  without  including  the  center  section. 

In  the  first  area,  the  stresses  as  found  in  Table  IX  show  fairly 
close  agreement.   In  the  panels  where  the  best  stress  agreement  occurs, 
such  as  panel  38,  the  strain  readings  on  both  sides  of  the  skin  are 
nearly  equal.   The  strain  readings  of  back-to-back  gages  were  averaged 
and  used  to  find  the  experimental  stresses  at  each  station.   Where  the 
corresponding  interior  and  exterior  gages  showed  considerable  variance 
in  their  strain  readings,  the  experimental  and  theoretical  stresses  also 
had  greater  variance.   The  difference  in  strain  readings  is  due  to 
differential  bending  in  the  cover  skin  and  possible  torsional  resistance 
of  the  pivot  rib  at  the  rear  fitting,  both  of  which  were  neglected  in  the 
theoretical  analysis.   The  differential  bending  was  neglected  because 
the  model  used  for  the  flexibility  coefficient  of  the  cover  skin  was  a 
flat  panel  which  carries  only  shear.   The  torsional  resistance  of  the 
pivot  rib  was  neglected  because  the  effect  of  an  unknown  amount  of  slack 
movement  about  the  wing  hinge  could  not  be  incorporated  into  the 
theoretical  analysis.   The  strain  gages  themselves  are  believed  to  be 
accurate  to  within  5  micro-inches  of  the  true  strain  values.   Strain 
readings  found  by  Holgren  and  Comfort  in  1963  were  very  close  to  those 
found  in  this  study.   The  substitution  of  aluminum  pads  between  the 
jig  and  the  loading  floor  for  the  original  plywood  pads,  as  recommended 
in  Reference  2,  does  not  seem  to  have  made  much  difference  in  the  strain 
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values.   The  original  gages  were  carefully  installed  and  waterproofed 
which  seems  to  account  for  their  longevity.   Most  of  the  gages  have 
been  in  place  for  over  four  years.   The  normal  expected  life-span  of 
this  type  gage  in  industrial  use  is  one  or  two  years.   The  first  factor 
mentioned,  that  of  differential  bending  and  torsional  resistance  of  the 
pivot  rib,  seems  most  likely  to  be  the  cause  of  the  difference  between 
theoretical  and  experimental  results.   The  differences  do  not  seem  to 
be  too  great  since  some  of  the  references  mention  that  errors  of  100% 
are  sometimes  found  in  the  case  of  thin  skinned  structures.   More  con- 
clusive comparisons  could  be  made  if  the  torque  loading  were  increased 
and  larger  strain  readings  were  available.   However,  it  was  decided  that 
a  significant  increase  in  the  torque  loading  beyond  the  336,000  in.  lbs. 
used  in  this  study  might  permanently  warp  the  wing  and  render  it  useless 
for  further  study. 

The  second  area  of  comparison  is  between  the  theoretical  results 
found  in  this  study,  and  those  of  Reference  2.   Including  the  center 
section,  from  Y  =  0  to  Y  =  25.095  resulted  in  changing  25  flexibility 

WW 

influence  coefficients  and  did  not  change  the  equilibrium  equations.   In 

effect,  it  moved  the  wing  root  from  the  pivot  rib  at  Y  =  25.095  to  the 

w 

center  rib  at  Y  =0.   Most  of  the  coefficients  changed  by  3  or  4%,  and 
none  of  them  changed  by  more  than  1370;  yet  the  matrix  of  resultant  loads, 
j*AT/  ,  shows  some  considerable  changes.   In  comparing  (at|  of  Reference 
2  with  the  /AT]  matrix  found  in  this  study,  the  former  loads  will  be 
designated  as  aj  (Reference  2).   The  greatest  change  was  in  the  load 
on  the  rear  intermediate  beam,  R.I.B.   a.  (Reference  2)  =  1010.4;  and 
a,„  =  10,263.5.   This  load  is  applied  to  element  24  and  will  influence 
the  stresses  on  panels  23  and  25.   The  flexibility  coefficient  associated 
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with  this  force  has  a  much  smaller  change,  F,„,  ,_  (Ref.  2)  =  13.3832 
and  F,„  ,L   =    14.3406.   The  purpose  of  mentioning  these  changes  is  to 
show  that  the  cumulative  effect  of  minor  changes  in  25  out  of  129  non- 
zero flexibility  coefficients  may  be  great.   This  means  that  small  errors 
in  deciding  upon  the  true  area  and  length  of  the  elements  of  the  ideal- 
ized model  may  be  significant.   When  the  actual  structure  to  be  idealized 
has  irregular  cut-outs  and  joints,  and  reduction  in  skin,  web,  and 
flange  areas  occurs  in  steps  rather  than  linearly,  considerable  effort 
will  be  involved  in  just  constructing  an  idealized  model.   By  examining 
Table  X,  it  is  seen  that  inclusion  of  the  center  section  relieved  the 
stresses  at  the  leading  edge  of  the  wing  which  was  the  area  of  greatest 
discrepancy  in  Reference  2.   Both  normal  and  shear  stresses  are  much 
closer  to  the  experimental  values.   The  stresses  in  the  rear  panels  are 
also  improved,  but  it  is  felt  that  the  pins  at  point  B  and  C  in  Figure  2 
fit  snugly  enough  to  offer  some  torsional  resistance  of  the  pivot  rib  at 
the  rear  fitting,  which  was  not  accounted  for  in  the  theoretical  analysis. 
Tables  XI  and  XII  list  the  resultant  loads,  a.  through  a,Q  for  the  wing  with 
and  without  center  section  included.   The  agreement  between  theoretical 
and  experimental  stresses  throughout  the  wing  seems  to  bear  out  the  val- 
idity of  the  matrix-force  method, and  the  approach  as  described  in  this 
study  in  the  prediction  of  stresses  in  a  wing  prior  to  its  construction. 
For  greater  accuracy,  it  would  be  necessary  to  use  more  elements  which 
would  require  the  calculation  of  additional  equilibrium  equations  and  hence 
expansion  of  the  matrices. 
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CONCLUSIONS 

The  matrix-force  method  of  analysis  using  the  strain  energy  theory, 
as  proposed  in  Reference  1,  is  acceptable  for  a  thick-skinned  wing. 
It  can  be  expected  to  predict  the  stresses  in  a  complex,  highly  redun- 
dant structure,  loaded  in  torsion,  with  a  reasonable  degree  of  accuracy. 
For  increased  accuracy,  the  long  slender  panels  representing  the  cover 
skin  should  be  broken  up,  and  smaller  elements  used  which  would  reduce 
the  size  of  the  structural  grid  but  increase  the  number  of  unknowns. 
Reference  7  mentions  a  method  of  accomplishing  this  by  dividing  each 
panel  into  a  number  of  narrow  longitudinal  strips  which  could  be  treated 
as  rectangles  with  constant  shear  flow  in  each  one.   The  difference  in 
shear  forces  of  adjacent  strips  would  be  taken  by  assuming  concentrated 
stringer  areas  in  the  spanwise  direction.   This  method  would  still 
neglect  Poisson's  effect.   From  the  results  of  this  study,  it  appears 
that  a  better  panel  would  be  one  whose  length  and  width  were  of  similar 
size. 

In  summation,  it  was  found  that  the  accuracy  of  the  matrix-force 
method  was  improved  considerably  (compared  with  the  results  found  in 
Reference  2)  by  the  inclusion  of  14  additional  members  which  resulted 
in  changes  to  25  flexibility  coefficients.   Generally  speaking,  it  must 
be  possible  to  attain  a  very  high  degree  of  accuracy  by  fundamentally 
the  same  method  as  has  been  used  in  this  study,  by  dividing  the  structure 
into  a  much  larger  number  of  elements.   This  involves  a  larger  number  of 
redundancies  and  much  more  numerical  work.   This  could  be  programmed  to 
be  accomplished  by  the  computer.   An  experienced  stress  analyst,  who 
knows  the  stress  distribution  of  most  of  the  elements  of  the  structure 
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on  which  he  is  working,  would  be  able  to  concentrate  the  redundancies 
to  the  unknown  points,  thus  reducing  the  computational  work  considerably. 
The  experimental  analysis  would  be  improved  by  the  addition  of  more  strain 
gages  on  the  inside  of  the  wind  cover  skins. 

Because  of  the  number  of  man  -  hours  involved  in  removing  and 
replacing  the  cover  skins,  however,  it  would  not  be  practical  to  install 
additional  gages. 
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TABLE  III 
THE  EDGE  SHEAR  FORCES  ON  THE  PANELS 


Piece 

wi 

wc 

lb 

Id 

ab/ai  = 

ac/ai= 

ad/ai= 

No. 

in. 

in. 

in. 

in. 

xb/wc 

wi/wc, 

Vwc 

1 

5.64 

6.82 

14.3408 

14.3287 

2.1027 

.8269 

2.1009 

2 

6.82 

7.44 

13.5077 

13.5041 

1.8155 

.9166 

1.8150 

3 

7.44 

7.66 

13.0307 

13.0303 

1.7011 

.9712 

1.7010 

4 

7.66 

7.54 

12.5756 

12.57  54 

1.6678 

1.0159 

1.6678 

5 

7.54 

7.06 

12.1370 

12.1296 

1.7191 

1.0679 

1.7180 

6 

7.06 

6.30 

11.8046 

11.7984 

1.8737 

1.1206 

1.87  27 

7 

5.64 

7.24 

79.3021 

79.3021 

10.9533 

.7790 

10.9533 

8 

6.82 

8.94 

74.9206 

74.9206 

8.3803 

.7628 

8.3803 

9 

7.44 

9.66 

71.3846 

71.3846 

7.3897 

.7701 

7.3897 

10 

7.66 

9.70 

68.5890 

68.5890 

7.0710 

.7896 

7.0710 

11 

7.54 

9.20 

66.4567 

66.4567 

7.2235 

.8195 

7.2235 

12 

7.06 

8.32 

64.9038 

64.9038 

7.8009 

.8485 

7.8009 

13 

6.30 

7.20 

63.8339 

63.8339 

8.8658 

.8750 

8.8658 

15 

14.3287 

20.8705 

79.3021 

74.9206 

3.7997 

.6865 

3.5897 

17 

13.5041 

19.0991 

74.9206 

71.3846 

3.9227 

.7070 

3.7375 

19 

13.0303 

17.7  87  9 

71.3846 

68.5890 

4.0130 

.7325 

3.8559 

21 

12.5754 

16.5114 

68.5890 

66.4567 

4.1540 

.7616 

4.0249 

23 

12.1296 

15.2393 

66.4567 

64.9038 

4.3608 

.7959 

4.2589 

25 

11.7984 

14.1418 

64.9038 

63.8339 

4.5894 

.8342 

4.5138 

28 

11.8046 

14.1748 

64.9038 

63.8339 

4.5788 

.8327 

4.5033 

30 

12.1320 

15.2456 

66.4567 

64.9038 

4.3590 

.7957 

4.2571 

32 

12.5756 

16.5132 

68.5890 

66.4567 

4.1535 

.7615 

4.0244 

34 

13.0307 

17.7879 

71.3846 

68.5890 

4.0130 

.7325 

3.8559 

36 

13.5077 

19.1022 

74.9206 

71.3846 

3.9220 

.7071 

3.7369 

38 

14.3408 

20.8878 

79.3021 

74.9206 

3.7965 

.6865 

3.5868 
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TABLE  V 
MATRIX  B     INPUT  FORM 


11     1.0  134 

3  3  -1.0  335 

5  5  -1.0  537 

7  7  -1.?  8  1 

8  4  96.1099  8  5 
8  8  7.4179  8  9 
812  1.8124  813 
816  3.0187  817 
820  2.5220  821 
824  3.3099  825 
828  3.2013  829 
832  2.8758  833 
836  -45.6107  837 
840  -3.5  384  841 
844  -2.5220  84S 
848  -3.3090  84^ 
945  -1.0  1022 

1123     1.0  1136 

1237     7.2236  1248 

1349  -1.0  14  7 

1538  -7.8009  1543 

1642     1  .0  1730 

1829  -1.0  1835 

1934  -8.3804  1939 

2010  -0.99999  2011 

2014     0.99915  2015 

2018     0. 99980  201O 

2022  -0.60391  2023 

2026  -0.45373  2027 

2030     0.55916  2031 

2039  -0.65009  2  040 

2043  -0.48235  2044 

2047     0.55916  2048 

2121  -0.75927  2122 

2125  -0.87645  2126 

2129     0.79688  2130 

2133     0.89114  2139 

2142  -0.85597  2143 

2146     0.79688  2147 

22  1  56.885  22  2 

22  5  56.885  22  6 

22  9  1.12519  2  210 
2213  -2.34025  2214 
2217  0.27134  2218 
2221  3.39403  2222 
2225  3.64622  2226 
2229  3.84895  2230 
2233  3.2J813  2239 
2242  -3.93750  2243 
2246  -3.R4895  2247 

23  8  56.8554  23  9 
2312  56.8702  2313 
2316  -56.8827  2317 
232^  0.0  2321 
2324  -63.5718  2325 
2328  15.8462  2329 
2332  78.4498  2333 
2341  -47.9'j19  2342 
2345  15.3462  2346 
2349  78.4498  2412 
27  9  7.5385  28  8 
29  3  1.0  29  4 
29  7  1.0  29  3 
2911  -0.D0844  2912 
2^15  -0.02295  2916 
2919  0.03219  2920 
2923  0.01487  2924 
2927  0.01009  2929 
2931  0.01249  2932 
2935  -0.77119  2936 
2939  -0. 01415  2940 
2943  -0.00971  2944 
2947  -0.01487  294* 
3018  -1.0  3m 
3216  -1.0  331? 


O.J 

220 

1.0 

244 

-1.0      1 

1.0 

4  4 

-1.0 

436 

1.0      1 

1.0 

6  6 

-1.0 

63  8 

1.0      1 

55.2466 

8  2 

69.5754 

8  3 

8  3.0796   1 

108.6853 

8  6 

120.8150 

8  7 

132.6135   1 

5.9191 

810 

4.2884 

811 

3.0187   1 

0.5448 

814 

0.5448 

815 

1.8124   1 

4.2884 

818 

5.9191 

819 

7.4172   1 

3.2013 

822 

3.5384 

823 

3.5707   1 

2.8758 

826 

2.3642 

827 

2.5220   1 

3.5384 

8  30 

3.5707 

831 

3.3099   1 

2.3642 

834 

-15.9213 

835 

-30.7843   1 

-60.8687 

838 

-75.9532 

839 

-3.2013   1 

-3.5707 

842 

-3.3099 

843 

-2.8758   1 

-3.2013 

846 

-3.5384 

847 

-3.5707   1 

-2.8758 

921 

1.0 

934 

8.3804   1 

1.0 

1035 

7.3897 

1046 

-1.0      1 

7.0710 

1147 

-1.0 

1224 

1.0      1 

-1.0 

1325 

1.0 

1338 

7.8009   1 

0.0 

14  8 

7.^563 

1532 

-1.0      1 

1.0 

1631 

-1.0 

1637 

-7.2236   1 

-1.0 

1736 

-7.0710 

1741 

1.0      1 

-7.3897 

1840 

1.0 

1928 

-1.0      1 

l.n 

20  8 

-0.99948 

20  9 

-0.99980  1 

-0.99996 

2012 

-0.99974 

2013 

-0.99915  1 

0.99974 

2016 

0.99996 

2017 

0.99999  1 

0.99948 

2020 

-0.69666 

2021 

-0.65009  1 

-0.55916 

2024 

-0.51788 

2025 

-0.48235  1 

0.69666 

2028 

0.6  5009 

2029 

0.60391  1 

0.51788 

2032 

0.48235 

2033 

0.45373  1 

-0.60391 

2041 

-0.55916 

2042 

-0.51788  1 

0.69666 

2045 

0.65009 

2046 

0.60391  1 

0.51788 

2049 

0.48235 

2120 

-0.71732  1 

-0.79688 

2123 

-0.82936 

2124 

-0.85597  1 

-0.89114 

2127 

J. 71732 

2128 

0.75927  1 

0.82936 

2131 

0.85597 

2132 

0.87645  1 

-0.75927 

2140 

-0.79688 

2141 

-0.82936  1 

-0.87645 

2144 

0.71732 

2145 

0.75927  1 

0.82936 

2148 

0.85597 

2149 

0.87645  1 

56.885 

22  3 

56.885 

22  4 

56.885    1 

56.885 

22  7 

56.885 

22  8 

1.83113  1 

0.27134 

2211 

-0.4  8011 

2212 

-1.30551  1 

-2.34025 

2215 

-1.30551 

2216 

-0.48011  1 

1.12519 

2219 

1.83113 

2220 

2.59681  1 

3.84895 

2223 

4.02233 

2224 

3.9375J  1 

3.20813 

2227 

2.59681 

2228 

3.39403  1 

4.02233 

2231 

3.93750 

2232 

3.64622  1 

-3.39403 

22  40 

-3.84895 

2241 

-4.02233  1 

-3.64622 

2244 

-2.59681 

2245 

-3.39403  1 

-4.^2233 

2248 

-3.93750 

2249 

-3.64622  1 

56.8736 

2  310 

56.8844 

2311 

56.8827   1 

56.8367 

2314 

-56.8367 

2315 

-56.8702   1 

-56.8844 

2318 

-56.8736 

2319 

-56.8554   1 

-15.8462 

2322 

-31.8510 

2323 

-47.9019   1 

-78.4498 

2326 

-92.3667 

2327 

0.0      1 

31.8513 

2330 

47.9019 

2331 

63.5718   1 

92.3667 

2339 

-15.8462 

2340 

-31.8510   1 

-63.5718 

2343 

-78.4498 

2344 

0.0     1 

31.8510 

2347 

47.9019 

2348 

63.5718   1 

6.8182 

2511 

7.4397 

2610 

7.6599   1 

1.0 

29  1 

1.0 

29  2 

1.0      1 

1.  J 

29  5 

1.0 

29  6 

1.0      1 

0.03219 

29  9 

0.01978 

2910 

0. 00477  1 

-0.02295 

2913 

-0. 04114 

2914 

-0.04114  1 

-0.00844 

2917 

0.00477 

2918 

0.01978  1 

0.01009 

2921 

0.01415 

2922 

0.01555  1 

0.  )1249 

2925 

0.00971 

2926 

0.00705  1 

0.01415 

2929 

0.01555 

2930 

0.01487  1 

0.00971 

2933 

0.00705 

2934 

-0.76286  1 

-0.78969 

2937 

-0.81957 

2938 

-0.84856  1 

-0. J1555 

2941 

-0.01487 

2942 

-0.01249  1 

-0.01009 

2945 

-0.^1415 

2946 

-0.01555  1 

-0. 31249 

2949 

-0.00971 

30  9 

1.0      1 

1.0 

3117 

-1.0 

3211 

1.0      1 

1.0 

3315 

-1.0 

2 
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TABLE'  VII 

MATRIX'     F' 

Input  Form 

1 

1500.00 

2 

252,80 

3 

714.00 

4 

528.00 

5 

109.00 

6 

-210.50 

7 

-2782.00 

8 

-4523.56 

9 

-5212.17 

10 

-5568.74 

11 

-5395.86 

12 

-4593.58 

13 

-4654.94 

14 

-4654.94 

15 

-4593.58 

16 

-5395.86 

17 

-5568.74 

18 

-5212.17 

19 

-4523.56 

20 

2390.73 

21 

4825.58 

22 

5395.02 

23 

6125.17 

24 

6229.25 

25 

5802.52 

26 

2738.27 

27 

2390.73 

28 

4825.58 

29 

5395.02 

30 

6125.17 

31 

6229.25 

32 

5802.52 

33 

2738.27 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

2 
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TABLE  VIII 

EXPERIMENTAL  RESULTS  AT 

Yw  =  34. 
w 

5 

PANEL 
NO. 

GAGE  NO. 

STRAIN 
READINGS 

AVERAGE 
STRAINS 

O"  max 

Tin  ax 

Inte- 
rior 

Exte- 
rior 

Inte- 
rior 

Exte- 
rior 

CT   min 

15 

486 

61 

61 

1094.02 

485 

118 

118 

1008.22 

484 

-50 

-50 

-922.41 

21 

27  9 

459 

17 

25 

21 

851.68 

280 

458 

95 

100 

97.5 

726.87 

281 

457 

04 

-14 

-5 

-602.06 

25 

367 

468 

16 

-43 

-13.5 

678.15 

368 

467 

-47 

-99 

-7  3 

627.44 

369 

466 

29 

11 

20 

-576.74 

28 

28 

474 

-19 

34 

7.5 

960.47 

29 

473 

137 

86 

111.5 

827.86 

30 

47  2 

27 

-8 

9.5 

-695.25 

30 

89 

477 

9 

9 

9 

457.66 

90 

47  6 

55 

63 

59 

512.27 

91 

475 

-12 

-20 

-16 

-566.87 

34 

19 

480 

-10 

-21 

-15.5 

599.43 

20 

479 

76 

65 

70.5 

568.22 

21 

47  8 

29 

10 

19.5 

-537.02 

38 

13 

483 

-80 

-7  6 

-7  8 

495.37 

14 

482 

61 

62 

61.5 

803.49 

15 

481 

35 

42 

38.5 

-1111.62 
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Comj 

T 

^BLE  IX 

at  Station  Yw  =  34 

>arison  of  Results 

.5.   Load 

ing 

336, 

000  in.  lbs.  of  torque.   Principal  axis  is  measured 

from  load  reference  axis,  the  Center  Intermediate 

Beam, 

which  has  33°  58'  sweep  back. 

PANEL 
NO. 

MAX.   or'  MIN. 

normal;  stress 

MAXIMUM  SHEAR 
STRESS 

PRINCIPAL  AXIS 

Experi- 

Theore- 

Experi- 

Theore- 

Experi- 

Theore- 

mental 

tical 

mental 

tical 

mental 

tical 

15 

1094.02 

1605.47 

1008.22 

1333.88 

31°  52' 

30°  12' 

17 

1530.10 

1092.17 

28°  16' 

19 

1047.07 

895.67 

38°  32' 

21 

851.68 

678.85 

726.87 

845.91 

40°  52' 

40°  43' 

23 

-7  92.63 

819.44 

48°  05' 

25 

-602.06 

-970.55 

627.44 

1033.12 

38°  48' 

49°  20' 

REAR  BEAM 

28 

960.47 

767.42 

827.86 

870.20 

44°  44' 

47°  36' 

30 

457.66 

626.57 

512.27 

556.54 

39°  20' 

45°  24' 

32 

813.00 

669.19 

40°  12' 

34 

599.43 

851.95 

568.22 

761.82 

37°  26' 

39°  55' 

36 

-7  67  .  60 

901.15 

35°  49' 

38 

-1111.62 

-1052.83 

803.49 

1160.03 

27°  11' 

33°  56' 
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TABLE  X 

COMPARISON  OF  THEORETICAL  RESULTS  AT  STATION  Yw 

=  34.5 

(1) 

WITHOUT  INCLUDING  CENTER  SECTION, 

AS  IN  REF.  2. 

(2) 

INCLUDING  CENTER  SECTION 

MAXIMUM 

MINIMUM 

MAXIMUM 

PANEL 
NO. 

NORMAL  STRESS 

NORMAL  STRESS 

SHEAR 

STRESS 

(1) 

(2) 

(1)        (2) 

(1) 

(2) 

15 

2234 

1605 

-1280 

-1062 

1757 

1334 

17 

2025 

1530 

-651 

-654 

1338 

1092 

19 

1420 

1047 

-630 

-744 

1025 

896 

21 

1009 

67  9 

-67  9 

-1013 

844 

846 

23 

827 

846 

-619 

-793 

723 

819 

25 

994 

1095 

-418 

-970 

706 

1033 

28 

813 

97  3 

-373 

-7  67 

593 

870 

30 

714 

626 

-478 

-486 

596 

556 

32 

835 

813 

-505 

-525 

670 

669 

34 

1068 

852 

-47  8 

-67  2 

773 

762 

36 

1575 

1034 

-575 

-7  68 

107  5 

901 

38 

1982 

1267 

-1004 

i 

-1053 

1493 

1160 
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TABLE  XI 

MATRIX  OF  ALL  INTERNAL  LOADS 
{AT}-  [S]  {P} 
FROM  IDEALIZED  MODEL  WITH  CENTER  SECTION  INCLUDED 


1  174.64054260 

2  1325.450S2773 

3  462.C38330C8 

4  474.22583008 

5  675.57104492 

6  1:66.64038086 

7  -597.00024414 

8  -3451.01147461 

9  -34<38.  75415039 
1C  -4256.81250000 

11  -5356..  39062500 

12  -5181.C39C6250 
1?  -4892.48437500 

14  -73  50.367187  50 

15  -7291.81250000 

16  -5  829. 66 796 « 75 

17  -5258.5  1953125 

18  -4711.71093750 

19  -4013.84619141 

20  -  34:7. 25048828 

21  397;. 23242188 

22  2286.73657227 

23  2093.25366211 

24  -1755.84497370 
2  5  -5  79.25122^70 


26 
27 

28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
4  9 


266.98681641 

1781.35253906 

3944.2023925^ 

-8407.73046875 

-7190.99218750 

-5754.24218750 

-4264.46093750 

-2847. 58496094 

650.^9536133 

716.78784180 

392.5336^141 

-136.  54725647 

1679.07592773 

6767.  17968750 

6785.34375000 

5224.67187500 

985.44262695 

10263. 53515625 

3625.59106445 

61 34. 81640625 

510  3.66015  625 

3034. 79296875 

4810. 60546875 

6030.46484375 
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TABLE  HI 

MATRIX  OF  ALL  INTERNAL  LOADS 
{AT}«[S]{P} 

FROM  IDEALIZED  MODEL  WITH  CENTER  SECTION  NOT  INCLUDED 


1 
2 
3 
4 
5 

6 
7 
8 
9 
1C 
11 
12 
13 
14 
15 
16 
17 
18 
IS 
2C 
21 
22 
23 
2  4 
25 


-92.323eQ676 

1*92.31518555 

-6G. 23043823 

-222.62962341 

-253.n«P3483c 

-94.43901C62 

442.254636  67 

-2286.37548828 

-3785.3C151367 

-4119.  C 9765625 

-4763.9843750C 

-4985.  19921875 

-4969.66796675 

-4492.2734375C 

-435^.52343750 

-4  8  86,99609  375 

-5C1  ?.?195312c 

-  4  7  6  5  .  9  1  C  1  5  6  2  5 
-4219.2539C625 
-3867.  ^4985352 

5143.6R75QCCC 

2CC. 63256836 

-379.16625977 

-  1720.40966797 
-1354.33C56641 


26 
27 

28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
4C 
41 
42 
43 
44 
45 
46 
47 
48 
49 


-11CC.  95019531 

-939.57446289 

-859.61694336 

-3524.74560547 

-3361.27685  54  7 

-3227.84252C3C 

-2977.405C293O 

-2714.C4394531 

767.79956C55 

904.96679687 

592.9165C391 

-17.30769348 

-720.509C332C 

8591.921875CC 

8726.  140fc25CC 

7870.3671675C 

4969.57031250 

1C1C  .43847656 

546C.7343750C 

8681. 625 COCCC 

R357.C3125CCC 

798C.82C3125C 

4503.7617167? 

769.419A335<; 
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START 

CLEAR   ARRAYS 

FORM     A  ,   B     ,      UC 

1 

PATT      P-ATTQC      T 

RETURN 

APPENDIX  A 
PROGRAM  ENERGY 


YES  — 


PRINT  MATRIX 
SINGULAR 


INVERT      A     ,    T^ITE     AI 


CALL  MATMPY 
CJ  =  A"'B 


FORM     G     FRCM     CJ     ^Y 
ADDING  TRIVIAI   RELATIONS 
WRITE     G 


TRANSPOSE      G 


GAIT  MATMPY 
H  =  GrF  G  ,  WRITE  H 


FORM  P 


PARTITION  H 


FORM  H21 


KD 


FORM     '^22 
WRITE     H22 


CAI.L  GAUSS  3 


RETURN 


YESH  PRINT  MATRIX 
SINGUI AR 


'  '  2  T  "?°  ,  I  RITE  >P-?I 


3  )  (U 


CALL  MATMPY 
ENP  =  H22I  x  H21 


EN=  -ENP 


UTIN  —  I/EN 


CALL  MATMPY 
S  =  G  X  UNITN 
WRITE  S 


CALL  MATMPY 
AR  =  EN  X  P 

WRITE  AR 


CALL  MATMPY 
AT  =  S  X  P 

WRITE  AT 


FORM     C 
C-l/E  (  H11+H12  XEN   ) 
WRITE     C 


END 
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//ENERGY  JOB   ASZ 5,0 1FP , 5, 5  ,  MESSER SCHM IDT ,D.B .  ,MSGLEVEL  1 
//     FXEC  FORTCLG 
//FURT.SYSIN  DD  • 

C      COEFFICIENT  MATRICES  At  B.  AND  UC  ARE  BASIC  STRUCTURAL  DATA 
ODIMENSION  A(50, 50). AI  (50,50), B  (50,50), CJ  (50, 50).UC (50,50),  x 
1G(50,50).H(50,  5O)tGT(5O,5O),H22(5O,5O),H12(5O,50)  , UCG  (50,50), 
2H?  1(50,  50)  ,H22  1(50,  51),  EN  (16,33 )  ,  AR  (16)  ,  UNI  TN(50,  50)  ,S  (50,50)  , 
3AT(50),H12N(5O,SO)  ,C  (50,  50)  .LIST  (50),P(3  5)  ,ENP(16,33  ) 
OEQUIVALENCE (A, C J, UCG, H22 , H21 , S) , (B,G,H12,C),(UC,H), 
HAI.GT.H2  2I  ,UNITN,H12N) 
C      CLEAR  ARRAYS  TO  RECEIVE  INPUT 
DO  10  1=1,50 
00  10  J-l.50 
A(I  ,J-)  =  0.0 
B(I ,J  )  -  0.0 
UC(I, J)  =  0.0 
10  A I  (I,  J  )  =  0.0 
C      READ  NON  REDUNDANT  MATRIX  A 

31  CALL  INPUT(50,A) 
GO  TO  (31, 32), NEXT 

C      READ  LOAD  AND  REDUNDANT  MATRIX  B 

32  CALL  INPUT(5  J, B) 
GO  TO  (32, 33), NEXT 

C      FORM  ENEKGY  COEFFICIENT  MATRIX  UC 

33  CALL  INPUT(50,UC) 
GO  TO  (33,35), NEXT 

35  DO  36  I  =  1,49 
DO  36  J  =  1,49 

36  UC(J,I)  =  UC(T ,J) 

C      INVFRT  THE  A  MATRIX 

CALL  GAUSS3(33,1.F-U  ,A,AI,KER) 

GO  T0(365 ,364),KER 
364  PRINT  151 

151  F0RM4T(24H     MATRIX  A  IS  SINGULAR) 
36  5  PRINT  2  00 

37  PRINT  38 

38QF0RMAT(64H  A  MATRIX  INVERSE  BY  GAUSS3,  9  COLUMNS,  IE  1'9,  10  1 8 
1PER  PAGE// 
PRINT  100,  (f  AH  T,  J),  J- It  9)  ,  1-1,33) 
PRINT  20) 

PRINT  100, ((AI (I, J), J  =  10,  18), I  =  1,33) 
PRINT  200 

PRINT  1 00 i ((AI (I, J), J =19, 2 7), 1=1,33) 
PRINT  200  ,. 

PRINT  300,C(AI ( I, J), J=2  8,33), 1=1,33) 
DC  FORMAT(/9F10.6  ) 
150  FORMAT(/11F10.6) 
200  FDRMAT(lHl) 
C      POST  MULTIPLY  AI  WITH  B  MATRIX 

CALL  MATMPY(33,33,49,AI ,B,CJ) 
C      FORM  G  MATRIX  BY  ADDITION  OF  TRIVIAL  RELATIONS  TO  CJ 
DO  40  1=34,49 

40  C  J  (  I  ,  I  )  =  1  .  ■) 
DO  41  1=1 ,49 
DO  41  J=l,49 

41  G(I,  J)  =  CJClt  J) 
PRINT  200 

42  PRINT  43 

43  F0RMAT(50H   G  MATRIX  , FORMED  FROM  CJ  WITH  TRIVIAL  RELATIONS//) 
PRINT  100,((G(I,J  ),J  =  1,9),I=1,49) 

PRINT  200 

PRINT  100,((G(I,J),J=10,18), 1=1,49) 

PRINT  200 

PRINT    100 tCCGCIt J)t J=19,?7), 1=1,49) 

PRINT  20 J 

PRINT  100, (CG  (I ,J), J  =  28,36) ,  1  =  1,49) 

PRINT  200 

PR  I  NT  100,  (CG(I ,J), J  =  37, 45) , 1  =  1,49  ) 

PRINT  200 

PRINT  4  00, ((G  (I , J) , J  =  46,4  9),  1  =  1,49) 
300  FORMAT  ( /  6F13.7  ) 
400  FORMAT  (/  4F13 .7  ) 
C      TRANSPOSE  G 

DO  44  1=1,49 

DO  44  J=l,49 

44  GT(I,J)  =  G( J, I) 

C      FORM  H  MATRIX  BY  TRIPLE  PRODUCT  (GT)(UC)(G) 
CALL  MATMPY  (49, 49,49, UC,G, UCG) 
CALL  MATMPY  (49, 49,49, GT, UCG, H) 
PRINT  200 
PRINT  45 
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45 


50d 

6j0 


510 

57 


581 


54 

521 
52? 


524 
32  5 

53 


5M 
540 


5  6 

Oq 

62 
63 

nor 


81 

8 
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FORMAT 
PRINT 
PR  INT 
PRINT 
PP  INT 
PRINT 
PRINT 
PP  INT 
PRINT 
PR  INT 
PRINT 
PRINT 
FORMAT 
FORMAT 
RFAD  L 
DO  510 
P  11)     = 
RE  AD (5 
FORMAT 
GO  TO 
PARTIT 
DO  5  4 
DO  54 
I  I  =  I 
JJ  =  J 
H22CI  I 
PRINT 
PR  INT 
FORMAT 
PR  INT 
PRINT 
I  N  V  E  R  T 
CALL  G 
GO  TOC 
PR  INT 
PP  INT 
PR  INT 
PRINT 
DO  5  3 
DO  5^ 
KK  =  I 
H21CKK 
FORM  R 
DO  540 
00  540 
SUM  =  0 
DO  541 
SUM  - 
ENPCI 
DO  56 
TO  56 
ON  (  I  ,  J 
FORM  y 
DO  94 
DO  44 
UMITN 
DO  62 
UNIT'M 
DO  63 
DU  63 
K  K  = 
UN  IT  nIC 
FORM  S 
CALL  IM- 
PRINT 
PP  INT 
FORMAT 
PRINT 
PP  INT 
PR INT1 
PR  PIT 
PR  INT  1 
R FOUND 

un  8^ 

SUM  = 
DO  3  1 
SUM  = 
AH  CI  ) 
DO  59 
LISTCI 
PR  INT 


IX  H  FORMED  BY  TRIPLE 
,49) 


,J)  ,  J=l,9), I- 
iJ)  , J-10,18)  , 
,J),J* 19,27), 
,J  )  ,J>28,36), 
,J  ),  J=37,45;t 

,J  ),J=46,49), 

) 

4) 
X       P 


P(I)) .NEXT 
F12.4,     II    ) 

NEXT 
TR  IX 


(44H    MATR 

50?, re  hc  i 

20? 

500, ff H(I 

200 

500, ff H(I 

20u 

500,  CCHd 

200 

50),  CC  Hfl 

200 

60J,  ff  H(  I 

(  /9F13.4 

( /     4F13. 
OAD     MATR I 

I    =     1,35 

0 .  0 
,58)      I,     ( 

(    12       , 
f 57,581) , 
ION        H    MA 

1-34,49 

J  =  34,  49 
-33 
-33 

, JJ ) =  H( 
200 

522,  f(H?2 
C/8F12.4 ) 
2V0 

522  t  ((H22 
H2?  MAT 
AUSS  3  ( 1 
525,524) , 
15  j 
20  i 

522,  (Y  H22 
522,  (f H2? 
I    =  34,49 

J-1,33 
-33 

,  J  ")  H  (  I 

EOUNOAMT 

1=1,16 

K=l,  33 
.  u 

J  -  1  ,  1 6 
SUM    +■    H22IM  ,J)»H21f  J,K  ) 
K)    =     SUM 

1=1,16 

J =1,33 
)  =  -1. 
ATRIX   UN 
I 


PRODUCT,  GTXUCXG//) 


-1,49) 
-1,49) 
=  1,49) 
=  1,49  ) 
-1,49  ) 


I, J) 

fit  J), J  -  1,8), I  -  1,16  ) 


f  I, J) , J  =  °,16),  I  =  1, 16  ) 

R  I  X 

6,1.E-10,  H22,  H22I.KFR) 

KER 


IM,J),J  =  !,8),I  =  1,16) 
I  (I  ,  J  ), J=9,16) ,  1-1,  16] 


,J) 

LOAD  COEFFICIENT  MATRIX,  EN, 


(I 


=  1,33 
=  1,33 

,  J)  =   1 


»  ENP(I, J ) 

ITN  BY  PLACING  EN   BELOW   UNIT  MATRIX 


(I 


=  1,33 
, I)  =  1  . 
1=34,44 
J-1,33 

-33 

I  ,  J )  =  FN 
MATRIX  . 

A  T  M  P  Y  (49 

20.' 

80? 
f43H  STR 

15",  (fS(I 

20j 

US(i  , 


50, 
20. 

50, 
ANT 


usu , 

LOAD 


=  1,16 


33 


I 

r 

J  =  1 

SUM  +  EM 
=   SUM 
1=1,49 

)  =  ]  M 

?   ">  ! 


1 


CKK,J) 

S  =  (G)CUNITN) 
,49,33,  G,  UNITN,  S ) 


ESS  COEFFICIENT  MATRIX 
,J), J=l,ll ) , I =  1,49) 

J)  , J=12,2?) ,1  =  1,49  ) 

J  ),J»23,33), 1=1,49  ) 
MATRIX,  AR 


=  (G)  (UNITN)  //  ) 


(   I, J)*Pf J) 


*2 


PRINT  60 
600FORMAT  (50H  REDUNDANT  LOAD  MATRIX.   AR    -   H22I     H21     P   // 
1   16H  AR(  1-16)=  Af  )) 

PRINT  700,(LIST(J> ,AR(J) ,J=1,16) 
700  FORMAT  (  I  13X,  I  2  ,  F20.  8  ) 

TOTAL  UNKNOWN  INTERNAL  LOADS.   AT-  (S)(P) 

DO  640  1=1,49 

SUM-  0.0 

DO  641  J=l,33 
641  SUM  =  SUM  +  S(I,J)«P(J) 
64^  AT( I)  =  SUM 

PRINT  200 

PRINT  900  ,  . 

9000FORMAT  (  45H   MATRIX  OF  ALL  INTERNAL  LOADS.   AT »   (S)(P)   // 
1  16H  AT  (I)    -  A  C   )    ) 

PRINT  1000,  (L  1ST  CI  ),ATM)  t  1  =  1,49) 
1000  FORMAT  ( /13X,I2»F20*8) 

DO  52   1=1,33 

DO  52   J=34,49 
52  H12CI.J-33)  =  H(If J  ) 

E  -  10300000.0  .   •  ■ 

INFLUENCE  COEFFICIENT  MATRIX.   C  -  (  Hll  +  (H12)  (EN))  /  E 

DO  90   I  -  1,33 

DO  90   K  =  1,33 

SUM  -  0.3 

DO    95       J    =    1,16 
95    SUM   =    SUM    t     H12(I , J)*EN(J,K) 
90    H12Nd,K  •)    =    SUM 

SCALE    -     10000. 

DO    70       1-1,33 

DO    70       J-1,33 
70    C  CI  ,  J  ^    =    CCH(ItJ)      +       H12N(I,J  ))/E)*SCALE 

PRINT    200 

PRINT    73 
733FORMAT    (45    H    MATRIX    OF    DFFLECTION    INFLUENCE    COEFFICIENTS       / 
135H       C         =C(H11    +CH12HN))/    E)*(SCALE)) 

PRINT  150,  ((C( I, J ) ,J=  1,11) , 1=1,33  ) 

PRINT  200 

PRINT  150,  ((C(I, J), J=12,?2), 1=1,33^ 

PRINT  200  , 

PRINT  153* CCCCIfJ) ,J=2  3,3  3) ,1=1,33) 

RETURN 

END  ,  . 

SUBROUTINE   GAUSS3(N,EP,AS,XS,KER) 

IMPLICIT  REAL»8(A-H),REAL*8(0-P  ) 

REAL»4  AS(5D,53),XSC50,5n)*EP 

DIMENSION  (A  5 0 , 50 )  , X (50 , 50  ) 

DO  1  I=1,N 

DO  1  J=1,N 

A(I,J)=AS  (I ,J"> 

1  X  ( I  ,  J  )  =■  0  .  0 
DO  2  K=1,N 

2  XCK,K)=1.C 
1^  DO  34  L=l ,N 

KP=0 

Z=0.0 

DO  12  K=L,N 

IF (Z-DABS(A(K,  L)))  11,12, 12 

11  Z=DABS  A  K,L 
K  P=  K 

12  CONTINUE 
IF(L-KP)  13,2>,  20 

13  DO  14  J=L,N 
Z=A(L,J ) 
A(L,J)»A(KP,J) 

1^  AfKP,J)=Z 

DO  15  J-l.N 

Z=X(L,J) 

X  (  L  ,  J  )  =  X  (  K  P  ,  J  •) 
15  X(KP,J)=Z 
20  IF(DABS(A(L,D)  -PP)5  ),50,  30 

30  IF(L-Nm,34,34 

31  LP1=L+1 

DO  36  K=LP1,N 

IF(  A(K,D)  32,36,32 

32  RATIO=-A(K,L)/A(L,L) 
DO  33  J=LP1,N 

33  ACK, J)= A(K, J)-RATIO»A(L, J) 
DO  35  J=1,N 

35  X(K,J)=X(K,  J")-RATIO*X(L,  J  ) 

36  CONTINUE 
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3  4 
40 


41 

42 
43 


51 


-  r 

i 


21 

ln 


1  •■' '. 
1 .  \ 


CUNT 
OP  A 

I  I-N 

un  4 

lf(  I 
I IP1 

DC  4 

s=s+ 

XCI  I 

KER  = 

UO  5 

Dfl  5 

XS(I 

kLTU 

KER^ 

RETU 

ENn 

SURR 

r  i  me 
on  l 

no  l 

SUM 
DC  ? 

SUM 

en, 

^FTU 

end 

SUBR 

I'IME 

ruui 

(  K(  S 

READ 
FUR*" 
U()  1 

IF  C 
NT 


INUE 
^  1=  1 
*1-I 

3  J  =  l 
1 

i-ri)4 
-ii  +  i 

2  K=I 
A  CI  I, 

,  J)  =( 

1 

1  1  =  1 

1  J=l 

t  J  )  = 

RN 

2 

RN 

HUT  in 
NSI  CM 
"  1=  1 
.  K=l 
=  ■>.  J 
?  J  = 
=  SUM 
K  )  = 


,N 


1,43,43 

I  P  1  ,  N 
K)»X(K,  i ) 

x(  ii, j)-s)/a( ii, m 

,N 
X(  I,  J) 


F  WATVPY  ( L ,  M,  N,  A,  R,  C  ) 

AO'>,5  0,  R(50,50),  C (50, 5^1 
,  L 
iN 


1  ,  M 

+  A  (  I  ,  J  )  • 
SUN 


er  j,k) 


CCN 
IF  ( 
IF  ( 
IF  ( 

IF  ( 
CO    i 


h.FTU 
PR  PI 

rruM 

STOP 
END 


RN 

OUT  IN 
MS  I  ON 
VALL'N 
)  ,  12) 
(  r> ,  2 
\T  (  4 
'i  1  =  1 
K(  I  ) 
INUE 
I  l.NF 
I2.Nl 
I3.NF 
I  4 .  N  E 

n  cii 

RN 

t  in 

AT      (' 


F     INPUT    (N, * ) 

A  (  N ,  N  )  ,  K  C  6  ) 
r.E(Kfl),Jl),fKC2),J2),(Kf3),J3),(K('4),Il)f 
,  (Kf*).  13) 

1)     I1,J1,A1,     I?,J2,A2,     I3,J3,A3,     I4,J4,A4,     NEXT 
(212,  F12.6),  I  1  ) 

•6 

.r,T.    p;  )  C-,0    ro    lr'j 

.  *■• )  a  c  1 1 ,  j  n  -  a  i 

.     0)    A(I2,J2)    =    A2 
.     •'  )    A  CI3,  J3  )    =    A3 
.     .'  )    A  (14,  JM  =  A4 
2),    NEXT 

, K  ,  II  ,  J 1  ,  A  1  ,  I  2  ,  J  2  ,  A  2  ,  I  3  ,  J  3  ,  A  3  ,  I  4  ,  J4  ,  A4 
K  =  ',FS/     IX,     4(-M2,E12.M) 


6)  i 


APPENDIX  B 
PROGRAM  RECROSE 


START 


1  =  0  ,    L  =  0 


1  =  1  +  1 


READ     N,E,U 


READ  GAGE  N0.,€;,^,e3 


CORRECT  FOR  GAGE  FACTOR 


A=  e^e3 


B-  €,  -«. 


c    2e£-k 


0=1.0  -  II 


H=1.0  -  U 


R  =  SQRT   (B*+C*) 


%U  -   (V2H)R 


°^v   -    5/2   (A/G+R/H) 


or       =  E/2   (A/G-R/H) 


PHIPRI=l/2   ARCTAN  2   (C/B) 


NO—  PHIPRI  =  PHIPRI  +  1.5707 


PHIPRI  =  PHIPRI/.017U5 


D  =   T    (C0S2#PHIPRI) 


cr  =  (cr      +(T    )/2+D 

/         v  /nax        mm' I 


CFm  fa     +<r      )/2-D 

3       V  m«T    mm  '/ 


7!  =   -  T      *SIN  (2.0  PHIPRI) 

o  max 


NO  —  PRINT  TITI.E 


YES 


PRINT  DATA 


Kl) YES-*  GO  TO  1 


NO 


END 


© 
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// 
// 
// 

c 
c 
c 
c 
c 


RECROSE  JUB   ASZ5,01FP 
EXEC  FORTCLG 

FORT.SYSIN  DO  * 

RECTANGULAR  STRAIN 
OF  PRINCIPAL  AXES, 
AND  STRESSES  ALONG 


MESSERSCHMinT,D.B.  ,MSGLEVEL  1 


RUSETTE  DATA  REDUCTION  FOR  THE  DET ERM I N AT  I CN 
PRINCIPAL  STRESSES,  MAXIMUM  SHEARING  STRESSES, 
THE  ROSETTE  AXES  1  AND  3. 
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R 
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w 
H 
9S 
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I 

IN 
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I 

F 
F 
F 
F 
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3F 
3 

3S 

3T 
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EAD 

-0 

'0 

*I+1 

EAD 

1-E1 

2-F2 

3  =  E3 

ORRE 

1-E1 

2  =  F2 

3=E3 

=  1.  i 

-1.9 

-El* 

=  E1- 
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